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Using the concept of crossing state and the formaUsm of second quantization, we propose a 
prescription for computing the density of arrivals of particles for multiparticle states, both in the 
free and the interacting case. The densities thus computed are positive, covariant in time for time 
independent hamiltonians, normalized to the total number of arrivals, and related to the flux. We 
investigate the behaviour of this prescriptions for bosons and fermions, finding boson enhancement 
and fermion depletion of arrivals. 
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I. INTRODUCTION 



A long standing issue in the theory and experiment of quantum mechanics has been that of measuring and formal- 
izing time observables. In the last two decades a substantial body of work has been produced clarifying theoretically 
and measuring experimentally quantities such as dwell times Q], tunneling times, or arrival times j^. In particular, 
many recent papers have challenged the classical work of AUcock, who denied the possibility of defining a quantum 
arrival-time concept ^, |^. In fact, these theoretical efforts and difficulties concerning arrival times have been 
essentially decoupled from the daily practice of many laboratories, where time-of-flight (TOF) methods are routinely 
used. One reason for such a divorce is that, in most cases, a classical analysis of the translational motion and the 
associated arrival-time distribution is sufficient. It is now the case, however, that the development of laser cooling 
techniques is bringing the quantum nature of the atomic dynamics to the fore, thus approaching the conditions for 
testing several proposed time-of- arrival (TOA) theoretical distributions in a regime that differs from the classical 
approximation. 

Yet another difficulty for a comparison and further interaction between experiment and theory is the absence, up to 
now, of a TOA theory for multiparticle systems. While the possibility of detecting individual atoms with nanosecond 
time resolution in specific TOF experiments is open , in the generic case the TOF spectra are produced by clouds of 
many particles that may interact with each other or /and with an external field. The aim of this paper xis to provide a 
quantum TOA theory which is applicable for the generic (one dimensional) multiparticle case using the formalism of 
second quantization (see for example 0), together with the crossing states introduced in and developed further in 
. We shall also portray several numerical examples to illustrate the phenomena of boson enhancement and fermion 
depletion of common arrivals. 



II. TIME OF ARRIVAL OF A SINGLE PARTICLE 



One of the major hindrances to the consideration of time observables in the framework of standard quantum 
mechanics was Pauli's theorem, which, simply put, states that no self-adjoint operator can exist that has canonical 
commutation relations with a self-adjoint bounded or semibounded hamiltonian, thus implying that the standard 
recipe associating self-adjoint operators to observables cannot work for time. 

Nonetheless, Aharonov and Bohm considered the motion of free particles as a clock to measure time, and introduced 
a time operator by symmetrizing the classical expression for the time when a particle, initially at the origin and with 
momentum p, passes point x, that is, t = mx/p. With a sign change this becomes the time of arrival at the origin of 
a free particle that, at time t = 0, is at position x with momentum p pO| If the arrival occurs at X rather than 
at the origin, then the corresponding "Aharonov-Bohm time-of-arrival operator" takes the form 



{X ^x)\+hx -x) 

p p 



(1) 



For all practical purposes, this expression fulfills all the properties one would expect of an operator associated with the 
observable quantity time-of-arrival for free particles on the line jl^ (for more details on this and the following topics, 
see 1^). It cannot be applied onto states with non vanishing zero momentum, which has been at times regarded as 
a drawback O, fil. In fact, this "difficulty" is perfectly physical, and mirrors the classical divergence of the time 
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of arrival when the particle's momentum tends to zero. It also explains how Pauli's theorem can be circumvented: 
Aharonov and Bohm's time operator Tj^^b is a maximally symmetric operator, therefore not self-adjoint. Other steps 
had to be taken before maximally symmetric operators and their concomitant POVMs (positive operator valued 
measures or generalized non-orthogonal resolutions of the identity) were understood physically, however. 

Although the faith in Pauli's theorem could have been slightly shaken by Aharonov and Bohm's proposal, the issue 
seemed to be settled after the important series of papers of Allcock l^] , which apparently put to rest all hope to 
obtain a sensible prescription for the quantum prediction of times of arrival. Even so, some adventurous souls kept 
on searching for alternative formulations within quantum mechanics. Kijowski in 1974 ||l5|] put forward a procedure 
to compute time-of-arrival probability densities for the free particle case in a purely axiomatic way (see also similar 
later work by Werner 

With the advent of a better understanding of positive operator valued measures (also known as generalized decom- 
positions of the identity or non-orthogonal measurements) |8[ [l^, |2^, the force of Pauli's argument was strongly 
diminished. In fact, it has been possible to show the relation between Aharonov and Bohm's time-of-arrival operator 
and Kijowski's distribution: they follow naturally one from each other |2^, |2^. Tab is the first operator moment 
of a POVM whose distribution function over a given state is the corresponding Kijowski distribution. Kijowski's 
distribution Il{t) of times of arrival at X for particles in a state \ip{0)) = \ip{t = 0)), in the free particle case, can be 
written as 

n(t,X)=^|(^(0)|t,a)p, (2) 

a=± 

where \t,a) are the two generalized eigenvectors of Tab with eigenvalue t, TAB\t,a) = t\t^ot) for a = ±. Note that 
the above mentioned domain "difficulty" of Tab does not apply to the bilinear functional 11 (i), which may be defined 
for arbitrary physical states regardless of their behaviour at p = . 

The fact that Tab is not a self-adjoint operator is clearly identified from the non-orthogonality of the complete 
basis a)}j a=±- ^^^^ eigenvectors are related to each other by means of the relation 

|t,a) = e*^(*"*')/'"'|i',a) (3) 
which assures the invariance of the distribution with respect to time traslations. In particular, 

K,a) =e^^*/''|i;„), (4) 

where we have used a special notation for the t = (generalized) eigenvectors or "crossing states", \vo) = \t = 0,a), 
where again a stands for either -I- or — (we will not be denoting explicitly the point of arrival X, which is part of the 
definition of these states, but it is always implied). 

In terms of these states we may rewrite Kijowski's distribution of times of arrival at the point X as 

n(t,x) = ^|(^(i)k.)p. (5) 

Q = ± 

This also suggests a rewritting of the distribution in terms of an operator for the density of arrivals at point X, 

u{t,x)^m)mx)m))- (6) 

Consider now the explicit form of the states \va) in momentum representation, 

(Pko) = (^)'^'©(«P)e~^''^/^ (7) 
V nm/ 

where 9(-) is Heaviside's unit step function. The correct correspondence of Il{t,X) with the classical case becomes 
now evident. If the non commutativity of position and momentum operators could be neglected, tt would correspond 
to the sum of the moduli of the fluxes that cross X from both sides. Also important is the fact that in a classical setting 
the corresponding dynamical variable provides the arrival distribution irrespective of the dynamics and interaction 
potentials. In other words, Eq. (|^) generalizes the free motion case in a natural and simple way for arbitrary 
interaction potentials, a task that could not be carried out using the original axiomatic procedure of Kijowski or by 
quantizing the classical time of arrival for each particular potential (the expressions are not analytically known in 
general and pose formidable ordering problems). 
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Underlying this rewriting of the time-of-arrival distribution a change of emphasis is to be found: whereas in Eq. 
the time-of-arrival distribution is obtained from the overlap of the initial wavefunction with the states associated 
with arrival at the instant t, be it from the left (a = +1) or the right (a — —1), in Eq. (H) it is obtained as the 
overlap of the evolved wavefunction with the constant states \va) that measure arrivals. The first point of view is, in 
a way, predictive: given the initial state of the particle, one can predict when the arrivals will occur. In the general 
case, with interacting potentials, this view may also be adopted with \t,a) given by Eq.(^, where the appropriate 
Hamiltonian is put in each case. From the second perspective, which could be termed "unconditional" , the arrival or 
otherwise of a particle at a; = is directly measured in physical space at every instant, using local definitions that are 
in no way conditioned by the different potentials in which the particles might be moving. This point of view, inspired 
by Wigner's formalization of the time-energy uncertainty relation , was advocated in |^ , where the properties 
of the crossing states \va) were examined, and Eq. (|^) was put forward as an expression of density of arrivals also 
for the case of interaction. In |9j we rewrote some other distributions that had been proposed in the literature for 
time-of-arrival distributions of particles in a potential (Q, later superseded by ||2^; see also |^) in terms of crossing 
states, and showed that those defined in Eq. (Q) were the only ones considered that led to classical correspondence 
with the properties expected of such distributions. 

Another particularly relevant aspect of the change of emphasis is that it helps to understand that Eq. (|^) need no 
longer be normalized to unity. In which case n(<, X) is to be understood as a density of arrivals of one particle: there 
might be a non zero probability for the particle never arriving at X, or, if the interacting potential were confining 
(such as the harmonic oscillator), recurrences would appear corresponding to many different arrivals. Notice that 
n(i, X) is a density of arrivals, not of first arrivals only. 



III. SECOND QUANTIZATION AND TIME OF ARRIVAL 



Even though TOF experiments with single atoms might be available in not too distant a future, we need to 
understand better how to predict time-of-arrival distributions for multiparticle systems. Most suited for such a 
purpose is the formalism of second quantization. One must first realize that the distribution of arrivals is a property 
of the same nature as the current density, or the kinetic enery, namely, it is obtained as the sum of "single particle" 
contributions, irrespective of the external or internal interactions affecting the A^-particle system. This is a key 
observation to discard outright, even for free motion, quantizations that would provide two-particle terms. 

Let Up and represent the annihilation and creation operators that respectively eliminate and create a plane wave 

of momentum p. Similarly, 'ip{x) and ip''{x) act on the vacuum disposing of and creating a particle at point x. The 
canonical commutation relations read 

[ap,al]^^ S{p~ q) , and ijj{x),-ip'' (y) ^ = S{x - y) , (8) 

where, as usual, [, ]± stands for the commutator in the case of bosons and for the anticommutator when fermions are 
involved. The position operator is written as 

dx X ^{x)^^^ (x) , 



and the inverse of the momentum operator as 

r+oQ 

P " = 



^ = dp- alup , 

J ~oo P 



from which the following form for a generalization of the time-of-arrival operator of Aharonov and Bohm might be 
inferred {X = 0): 



= ^-r^ / dxdpdqdr [((5(p — g) + (5(p — r)) aja^ + 2a^ajarap] . (9) 

In this expression one can recognize a one-particle component but also a two-particle one. As pointed out above, this 
leads us to discard this procedure, because of its unphysicality. 

At any rate. Tab is only valid for the free particle case, a further limitation of this route. The proper quantization 
procedure for the multiparticle case starts, as noted above, by recognizing the additive character of the time of arrival 
in terms of single particle contributions. 
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The basic trick is that for additive quantities taking the form of a sum of single particle operators, 

G = ?i+52+. •••+?JV, (10) 

each of which has matrix elements gji — {j\g\i) in a complete (single particle) basis, the multiparticle operator in 
second quantized form is given by the simple expression 

G = ^5jia]a, , (11) 

where Ui and aj are the i-th annihilation and j-th creation operators. That is, they connect states \i) and \j) 

respectively with the vaccum state {a.i\i) = |0) and \j) — a||0)). 

In the case of the arrival density operator we can directly apply this procedure in momentum representation. An 
even more compact expression is obtained by using the crossing states, to generate crossing operators, both annihilation 
and creation. Consider any (generalized) one particle state We can write the annihilation and creation operators 
associated with the state as 

/•oo poo 

/ dp {ip\p)ap ; (p^ = I dp{p\ip)al. 



On applying this procedure to the crossing states, we obtain the crossing operators 

'■ + °° r + °° / ^,r, \ 1/2 



MX) ^ / dp{vMap= ^P(^) 0(«P)e'^''/% ; (12) 



— OO J —oo 



^i(^) = J_ dpipMal^J^ '^pO e(ap)e-^'^/^at . (13) 

Let us now put together Eqs. ( p^ ) and ( p^ with Eq. ^ to write the arrival density operator n(X) for arrivals at X 
in second quantized form, 

n{X) - Yl ^iiX)v^iX) - 1" dpdq ^Q{pq)e^(^-P~>^/^ala, . (14) 
The left and right arrivals density operators 11+ (X) and n_(X) are similarly defined as 



We may also write from Eq. (14) the corresponding operator in Heisenberg picture, whose expectation value over the 
initial state will give us the density of arrivals at point X at instant t, 



u{t,x) = c/t(t)n(x)f/(t) = ;^^t(x,t)tj„(x,t) (15) 

= p dpdq^<d{pq)e^('i-^^'"'^Zl{t)Z,{t), (16) 

where a^(i) and ag(t) arc the time evolved creation and annihilation operators, with evolution operator U{t), i.e. 

at(t) = U^{t)ZlU{t) and similarly Zp{t) = U'^ {t)aj,U {t) . 

The density of arrivals at instant t at point X for a generic state |?/;) may thus be written as 

n(t,x; V) = {mm.x)\m) = mm{x)\m) ■ 

This expression agrees with Eq. (^ whenever \ip) \s a. one particle state. Even though it is not immediately apparent 
from expression (|l§) that we are obtaining positive semidefinite distributions, this is indeed the case by construction: 
n(t, X) is a positive operator because it is a sum of two terms of the form A. 

Furthermore, the one particle operator that is an extension of Aharonov and Bohm's time-of-arrival operator at 
position X for many particles, even in the interacting case, is straightforwardly written as 

fx - r dttn{t,X)^ f dtdpdq^Q{pq)e'^''^'P^^'^tal{t)aJt). 
J-oo J~oo hm 
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By construction this is simply a one particle operator, which coincides with Tab over states whose content is just one 
free particle. 

If the evolution of the system is governed by the free particle Hamiltonian it is easy to check that the integral over 
time of the arrival-density operator 11^'"^° {t, X) sums to the total particle number operator, 



/ + 00 I'+oc 
dt n'-^^^ {t,X)= dp dlap = N . 
-oo J — oo 



This is no longer the case whenever the evolution operator is not the free one; anyhow, we deduce from this expression 
that the arrival density is normalized to the total number of arrivals. Notice that in the interacting case the total 
number of arrivals need not coincide with the total particle number, it may be smaller or bigger. 

A particularly important property of the arrival-density operator is that the density of arrivals n(i, X; tp) over any 
state is covariant in time if the Hamiltonian is independent of time (as has been assumed all along). Even though 
the properties of covariance, positivity, and correct classical correspondence do not, by themselves, completely fix the 
density of times of arrival, they are minimal requirements, the lack of which would seriously impair any proposal. 

Even though in the presentation above we have restricted ourselves to pure states, there is no problem in extending 
our proposal to mixed states, as follows: 

n(i, X; p) = Tr (n(t, X)m) = Tr (fl{X)p{t) 

For the sake of completeness, let us note down the flux operator for many particles, in Schrodinger's picture, 

= ^ {iiHx)dx^^{x) - Idxi^Hx)] m)} = 

2m L L J J 

1 

^ 2h^J '^P'^'^ e^^^-P'^'"'^ {p + q)ala, . (17) 

or in Heisenberg's picture as 

^ ^ ^ ^ 1 f°° 

j{t,X) = uHt)jiX)Uit) - dpdqe^'^^-P^''^\p + q)dl{t)a,it), (18) 

(again assuming that the Hamiltonian is independent of time). Notice that the flux, defined in this standard manner, 
is a one-particle operator. 

A straightforward comparison of Eqs. (|l6|) and (|8|) reveals the differences and similarities between H and the flux. 
In the former a geometric mean of the momenta takes the place of the arithmetic mean in the latter. Moreover, H 
counts the case when p and q are both negative as a positive contribution to the arrival density, whereas the same 



case counts as a negative flux contribution in (18). This means that the quantity that tends classically to the flux is 
n+ — H_ rather than H itself. 

IV. FREE PARTICLES: BOSON ENHANCEMENT AND FERMION DEPLETION 

We have already made out several properties of the proposed arrival-density operator, namely positivity, covariance, 
one-particle status, classical limit, and normalization to total number of arrivals. There is an obvious missing element 
yet, in that we have not investigated so far whether the fermionic or bosonic character of the particles involved is 
somehow reflected in the properties of the distributions of times of arrival, as is to be expected. 

In fact, this distinction between fermions and bosons is already present in the proposed arrival-density distributions, 
as we will be showing in this section. In order to portray this new property it is enough to consider simply two-particle 
states, of generic form 

/oo I'OC 
dpidp2ip{pi,P2)\pi,P2) = / dpidp2ip{pi,P2)aldl\0) , 
-oo \ 2 J —oo 

both for bosons and fermions, where |0) is the vacuum state, and the normalization condition reads 

dpi dp2 



1p{pi,P2)±i}{p2,Pl) V(Pl,P2) = 1 
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where the upper sign corresponds to bosons and the lower one to fermions. 
Consider ^{pi,P2) given as 

V'±bl,P2) = - , , [Xa{pi)XbiP2) ± Xa{p2)Xb{pi)] , (19) 

which fulfills the normalization requirement if Xa and Xb are normalized one particle wavefunctions. Quite obviously, 
(XalXft) stands for J dp Xa{p)Xb{p) ■ In order to compare with the case of distinguishable particles, we shall also be 
using 

i^d{Pl,P2) = Xa{Pl)Xb{p2) ■ (20) 

Since the arrival-density operator Tl{t, X) is a one-particle operator, the density of arrivals over the state \'4'+), say, 
can be reorganized as 



Ni 



where N+ = ^2 (1 + |(xa|Xb) P), and 

U^,{t,X) = Y.{x^\vi{X,t)MX,t)\x,) 



a=± 



Over the fermionic state {ip-) the cross terms carry a negative sign in front. The evolved crossing states are given by 



Eqs. (|l^ ) and (13) on substituting cip and a], by ap{t) and (t), respectively. 

On the other hand, the evaluation of the expectation value of the evolved arrival density operator over the state 
IV'd), which computes the density of arrivals for two distinguishable particles in such a state, produces just the two 
diagonal terms, i.e. 

{MO)\^it, X)\M0)) - ^aa{t, X) + Ilbb{t, X) . 

It should be observed that these computations are general in that they hold true for the case of interacting particles 
as well, as long as the states have the form given above. These results indicate that fermions and bosons (antisymmetric 
and symmetric states) present cross terms in the density of arrivals completely analogous to those that in spatial density 
signal the statistics of the particles. In fact, the formalism of second quantization carries in itself the fermionic or 
bosonic character of the particles concerned, through the commutation relations. 

As a consistency check one may compute for the above states, iV'i) and \ipd), the corresponding reduced one particle 
density operators p^^\ j = 1,2 {p^^^ = p^^^ for IV'i)) and note that in all three cases n(t,X) =^^.n«(t,X), where 

n'-^\t,x)^TTj [n(t,x)p<^')^ 

in agreement with the one-particle character of the arrival-time distribution. 

The difference between bosons, fermions and distinguishable particles {\ipd)) is quite apparent in Fig. ^ The one 
particle states Xa and Xb are gaussians with a spatial separation between them (in atomic units, 3.5), while their 
width is 1 (a.u.). Correspondingly, there are two main arrival times (maxima of the arrival densities) for all three 
cases. Even so, the two sets of principal arrivals for bosons (symmetric state) are much closer together and much less 
differentiated than for distinguishable particles {\ipd)), which in turn present closer and less differentiated maxima 
when compared to the fermionic (antisymmetric) case. It should be noticed that in this situation of free motion the 
distributions are normalized to 2, as can be readily checked in this numerical simulation. 



V. INTERACTING PARTICLES 



Consider now a pair of interacting particles, be they distinguishable, bosonic or fermionic, moving in otherwise free 
space. The two-particle subspace of Fock space can be rewritten in center of mass and relative coordinates, and we 
shall consider for simplicity factorized states of the form 



■0(^1,^2) = X{P)Hp) : 
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FIG. 1: The solid line corresponds to distinguishable particles {\ipd}), the dots show the density of arrivals for bosons (symmetric 
state \ip+)), and the dashed lines that of fermions (antisymmetric state \Tp-)). The states are defined by Eqs. and pc|, with 
Xa and Xb gaussian states with minimum uncertainty at t — 0, their central positions being (Xal^lXa) ~ —3.5 and(xi)|S;|xf>) = 0. 
In both cases, their central positions in momentum space are at 3, and their spatial widths Ax = 1 (where Ax is the square 
root of the spatial variance). The point of observation is X = 3, and the mass m = 1. All magnitudes are expressed in atomic 
units. 



where P = pi + P2 is the center of mass momentum and p = {pi — P2)/2 the relative one. Under exchange of the 
particles P is unchanged, while p flips sign. So in order to ensure that the state is bosonic we are forced to use even 
functions (j)+{p) = 4'+{~p)j whereas the fermionic case demands odd functions 4'-{p) = —(/)_(— p). The total mass is 
2m, while the reduced mass ^ pertaining to the relative system is m/2. The normalization condition is translated 
into the requirement that 4>± and x(^) be normalized to unity. 

In what follows we shall assume that the center of mass function is gaussian with minimum uncertainty product at 
t = 0. As to the internal states, they will be evolving in a harmonic oscillator potential. We shall consider stationary 
and coherent internal states. 

Figs. H and || represent two different sets of cases concerning internal stationary states. The ground state and the 
even excited states are symmetric (bosonic), whereas the odd numbered excited states are antisymmetric (fermionic). 
The differences between Figs. ^ and ^ are due to the different ratios between internal energy and that of the center 
of mass motion. In Fig. |^ the internal oscillations are much slower than the center of mass motion, so the humps of 
the internal spatial wavcfunction appear, somewhat distorted for later times because of the spreading, in the arrival 
density. However those humps are smoothed over in Fig. js] due to the much slower center of mass motion relative 
to the internal motion. Correspondingly, the integral of the curves in Fig. ^ is very nearly 2, whereas there is a 
significant increase of this number in Fig. ^ for the excited states. The higher the excitation the broader the state is, 
spatially, thus leading to more crossings. 

We have also studied a case where the internal motion is time dependent. If it is fast enough with respect to the 
translational motion, a peak structure corresponding to several oscillations may be observed. Let us consider, at 
time t = 0, symmetric and antisymmetric combinations of coherent states, of the form (remember that the coherent 
state \z) is given by exp(— |2:p/2) ^^Q(z"/\/nI)|?^), where \n) is the n-th excited state of the harmonic oscillator 
hamiltonian) 

|0±) = i= [\z) ± \z)] , (21) 

where z is the complex conjugate of z. We shall take in particular z — i. In the relative motion space |z) is a minimum 
uncertainty product gaussian centered at the origin with average momentum (2/iti;?i)^/^ and spatial variance h/(2ujm). 
As time progresses it oscillates back and forth along the relative motion coordinate with period T = 2tt/uj. 

In Fig. H the translational motion is faster than the oscillations so we see just one peak for the symmetric case and 
two maxima for the antisymmetric (fermionic) case. Since there is hardly any component of negative momentum. 
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1 




t (a.u.) 

FIG. 2; Arrival density of two identical particles in an internal stationary state of the harmonic oscillator. The center of mass 
state is gaussian, initially centered on a; = and p — 4 with width Ax = 0.5. Circles correspond to the (internal) ground state, 
dashes to the first excited state, solid line to the second excited state, dot-dash to the third one. The internal frequency is 
u! — \/0.02 and the oscillation period is T ~ 44.4. The point of crossing is a; = 3. All magnitudes in atomic units. 



0.6 




-5 5 15 25 

t (a.u.) 

FIG. 3: As in figure |2| with the center of mass central momentum changed to 1, Ax = 1 and the internal frequency to lu = \f2. 
T ^ 4.4 

there is no distinction between flux and density of arrivals. 

In contrast to Fig. ^, the internal potential is much stronger in the situation depicted in Fig. ||. As opposed to the 
case of stationary internal states, the tighter binding produces here oscillations that can be clearly seen. 
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1.5 




-1 1 2 3 4 5 



t (a.u.) 

FIG. 4: Arrivals density (solid line for fermions and dashed line for bosons) and flux (triangles for fermions and circles for 
bosons) for the internal states defined in Eq. (^l|), with z — i and internal frequency uj — VO.02. The initial center of mass 
state is a gaussian with central position at a; = and central momentum p = 4. The width is Ax — 0.5. The point of arrival is 
X = 3. All magnitudes in atomic units. 
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-2 4 10 16 22 28 

t (a.u.) 

FIG. 5: Arrivals density (solid line for fermions and dashed line for bosons) for the internal states defined in Eq. (|2T|), with 
z = i and internal frequency lu — \/2. The initial center of mass state is a gaussian with central position at x = 0, central 
momentum p = 1, and spatial width Aa; = 1. The point of arrival is X = 3. All magnitudes in atomic units. 



VI. DISCUSSION 

In this work we have proposed a general method for computing densities of arrivals (and related arrival density 
operators) for multiparticle states, that fulfill a number of quite sensible demands: positivity, covariance (if the 
evolution is homogeneous in time), related to a one particle operator, normalized to the total number of particles 
in the free case, related to the flux, and consistent with the classical arrival density. The analysis of the density 
of arrivals (i.e., the one point function of the arrival density operator) shows consistency with the results one would 
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expect for bosons, fermions, and distinguishable particles. Numerical computations also show the behaviour expected, 
both for the free and the interacting case, and reveal a number of physical effects, hitherto unexplored. The proposed 
distribution is also applicable to the case of external interaction potentials as shown already for the single particle 
states. 

The fact that the arrival-density operator is a one particle operator implies that in fact no distinction is made for 
the one point function (the density of arrivals) between the bosonic case and the symmetric states of distinguishable 
particles. The full difference will be seen in two-point and higher order functions, such as the arrival - arrival 
correlation function. In fact, on computing this two point correlation function, one sees immediately that the one 
particle component of the two point operator behaves in the natural way and can be substracted from the arrival - 
arrival correlation function to give the correlation function of pairs of arrivals, {ip\ : II(t)II(0) : \tp), where : : stands 
for normal ordering. The full analysis of this object we will leave for future work. 

In this paper we have not been overly concerned with domain problems and the like. We know that the one particle 
time operator in the free case is only maximally symmetric and cannot be made self-adjoint, and we have no reason 
to expect that the interacting multiparticle case will be simpler in this respect. Nonetheless, this is not particulary 
relevant for our main interest, which lies in the computation of densities of arrival. Notice furthermore that the fact 
that Aharonov and Bohm's operator is not self-adjoint is really no hindrance to a full quantum mechanical analysis 
of its associated densities. 

There might be many other alternative prescriptions for times of arrival. In the present state of knowledge, we do 
not feel able to discard those outright. However, by pushing to the multiparticle case the definitions used for one 
particle, this analysis becomes more amenable to experimental test. 
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